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1.1  I Dl'.NTl  FI  (:AT10^■ 

V.  lu!:n  I.  Subniaiiiii  (iMF.dA  (lomput  t- 1'  I’ropiam  I’c  r 1 ormi  iicc  Specification,  ih’lino.s 
th.  liiiutional  m ipi  i ri 'r.icnt  lor  the  Submarine  OMI'KIA  domputi’r  I’roKram  whicli 
is  u.seil  b',-  the  AN/AKN-'»‘i  OMFKIA  Navij'.ation  Set.  The  Navipation  si't.  and  the 
()MKi:A  prop.r.im  topelher  cc'miiri.se  the  Suhmarine  OMFKIA  Kavijuition  Sy.stem.  The 
t .ipi  which  ibiiiie.s  the  computer  |)ro)',r.ain  i .s  entitled  AN/HKN-7  Navigation 
I’rogram. 

V.'lume  II,  Submarine  t'MliOA  (aimi)uler  Program  Design  Spi'c  i 1 i c a t i on  , .allocates 
the  lunction.il  re(|  u i remeiit  s of  Voluim'  1 to  the  computer  I'oul  i ne  .and  suh- 
prog.r.im  I eve  1 . 

This  volui’ie  de.stribes  the  suliprop,  r.ain  desip.nated  .as  Prop.ig.at  ion  Precl  i c- 1 i on  , 
which  h.as  the  .ihbre  vi  at  ion  PP  in  the  iirop.r.im  listing  (Volume  Xiil). 


\.l  PKOI’ACATION  PKKDf  CnoN  SrilPKOCPAM  TASKS 


1.2.1  Oyi  rvi  ew 

In  order  to  successfully  determine  the  current  submarine  position  from  the 
phase  information  received,  it  is  necessary  to  know  the  phase  velocity  of 
the  incoming  wave  along  the  path,  and  the  path  length.  To  obtain  a more 
accurate  estimate  of  .submarine  position,  we  will  investigate  a number  of 
effects  (magnetic,  day./night  (diurnal),  ground  conductivity,  and  variations 
in  radius  of  the  earth)  that  cause  position  to  deviate  from  that  calcu- 
lated using  constant  velocity  over  an  assumed  perfect  sphere. 

Operation  of  the  OMFGA  Navigation  System  is  based  on  the  measurement  of  the 
phase  of  several  transmitters  operating  In  the  10-14  kHz  electromagnetic 
spectrum.  It  has  t)een  established  that  the  eight  proposed  O.'-IEGA  trans- 
mitters will  cover  tfie  earth  with  signal  levels  adequate  to  permit  that 
phase  measurement.  Furthermore,  it  is  generally  accepted  that  signals  in 
this  spectrum  prop.agate  through  a waveguide  made  up  of  two  concentric 
spheres,  one  sphere  is  tfie  earth;  the  other  is  the  iijnosphere. 

Under  ideal  conditions  a simfile  phase  measurement  would  suffice  to  precisely 
locate  a receiving  station.  However,  die  walls  of  the  waveguide  are  not 
perfect  and  .ire  affected  by  several  parameters,  specifically  the  effect  of 
the  sun  on  the  ionosphere;  earth's  magnetic  field;  ground  conductivity,  and 
others.  Tliese  im.per  feet  ions  in  the  waveguide  walls  cause  changes  in  the 
propagation  of  tfie  (■•lectromagnet  ic  signal,  the  phase  velocity,  which  reduce 
the  accuracy  of  die  phase  measurement. 
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Tlie  phase  velocity  of  such  waves  in  a perlect  wavcipiuie  with  this  geometry 
has  been  found  to  depend  on  the  width  of  tlie  wavep,uide  (the  height  of  the 
ionospliere)  and  on  the  electrical  conductivity  of  the  surfaces.  The  elec- 
tric conductivity  of  the  earth's  surface  is  important,  as  well  as  that  of 
ionosphere.  • In  the  case  of  the  ionosphere,  the  problem  is  complicated 
theoretically  because  a continuous  wall  does  not  exist;  rather,  the  electron 
density  and  collision  trequency  vary  witli  height  in  some  manner  which  c.jn 
be  approximated  as  exponential.  The  effect  of  the  earth's  magnetic  field 
on  the  pliase  velocit;  must  also  be  considered.  This  has  been  analyzed 
theoretically  and  the  assymetry  between  propagation  from  cast  to  west  and 
from  west  to  east  has  been  observed. 

These  factors,  along  with  one  to  account  for  the  oblateness  of  the  earth, 
are  incorporated  into  a computer  program  designed  to  provide  incremental 
real-time  corrections  along  the  propagation  path. 


1 . d . J Tasks 

Tluri'  an-  lour  tasks  ol  the  Propagation  I’retliction  Subprogram.  They  are  re- 
1 1 Cl  t cd  i n I'i  g.urc  I . 

a)  1 111  t ial  i zal  ion:  (,'a  Icul  at  i on  of  those  environmental  parameters  wliich  are 
independent  of  ire(juencv  and  ni'ed  only  lie  calculated  once.  Included 
are  : 

• Tlie  central  patli  angle  lietween  station  and  OMI’.CA  receiver, 

• Till-  first  magnetic  parameter, 

• The  seasonal  Index,  and 

• Tlie  sun  vtu' t or  . 

b)  Integration  Over  the  Transmission  Path:  In  order  to  evaluate  the  effect 
ol  the  i.mosplure  on  the  transmitted  signal,  it  is  necessary  to  super- 
pose the  known  signal  path  on  the  model  of  the  world  contained  in  the 
suliprogram.  The  following  subtasks  are  frequency  or  patli  dependent; 

• Calculation  ol  the  diurnal  constants, 

• Auroral  compensation, 

• Karth  conductivity  Index, 

• Diurnal  function  and  average, 

• Correction  lactors  lor  propagation  modes,  and 

• Integration. 

c)  Generation  of  Output  Data:  As  a result  oi  the  incremental  integration 
the  subprogram  generates: 

• Averaging  constant;  correction  for  Incremental  integration, 

• Base  station  control  angle;  used  for  phase. 
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I PROPAGATION  fPfOICTI 
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iR.'CPr*.. /NT  or  patmj 


• start  inttopation 
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FIGURE  1 
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• I’ropagal  i on  CorrccLion,  9 . 

2 ^ 

• I’redicted  Variance  o , . 

PPi 

• Splu'i'ocdal  Correction,  9^. 

d)  Tracking  Kilter  Update:  Includes: 

• Predicted  Tracking;  Filter  measurement. 

• Pane  count  corrections. 
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SKCTION  2 

Al’Pl.lCABLK  DOCIIMKNTS 


a)  Submarine  OMKCA  Computer  Program  Periormanco  Spec i 1 LcaL i on  (Volume’  ( of 
llie  Submarine  OMKCA  Computer  Program  Specification) 

Applicable  Sections: 

5 . 1 Introdiic  t ion 

3.2  Functional  Description 

1.2.5  Detailed  System  Operations 

3.3  Detailed  Functional  Requirements 

3.1.10  Tracking  Filters 

3.3.11  Combinational  Kalman  Filter 

3.3.12  Propagation  Prediction 

b)  Submarine  OMF.CA  Computer  Program  Design  Specification  (Volume  11  of  the 
Submarine  OMI’ICA  Computer  Program  Specification). 

c)  NOKT  68-Ab,  NAP70  User's  Manual,  .July  1908. 

d)  NORT  h8-115A,  Detailed  Description  of  NDC-1070,  Computer  Instructions, 
Revision  A,  February  1970. 

e)  NORT  b9-H7A,  NDC-1070  Flow  Chart  Program,  User's  Manual. 
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SKCTTON  3 

rkquikkmi-;nts 


In  onltT  to  understand  the  sub -i)rogram  description  contained  in  the  follow- 
ing pages,  it  is  necessary  that  tlu'  reader  will  have  become  1 ami  1 i ar  with  the 
associateii  functional  requirements  found  in  Volume  I,  Performance  Specifica- 
tion, and  witli  t lu'  sub-program  allocation  found  in  Volunx"  11,  Design 
Spi-’c  i 1 i cat  ion  . 


i.l  DKTAII.IID  DKSCRIPTION 


5.1.1  Keferonce  Labels  to  Flow  Diagrams 

Tlie  code  used  to  reference  the  [wrticular  block  in  the  flow  diagraijis, 
Section  3.2,  is  as  follows:  The  first  number  is  that  page  number  found  in 
the  upper  ri  glit  corner  of  the  diagrams  . Tliis  will  be  followed  by  a slash 
sign  (/)  to  separate  the  page  number  from  the  block  designator.  Tlie  desig- 
nator will  eitlier  be  a mnemonic  label  (e.g.,  TKST  SYNC),  a local  label 
indicated  by  a dollar  sign  ($),  or  an  integer.  The  two  types  of  labels 
reference  the  particular  information  block,  on  the  given  page,  to  which  the 
lal)el  is  attached.  The  integer  number,  n,  means  that  the  referenced  block 
is  the  nf't  block  from  the  top  of  the  page;  p8/3  would  refer  to  page  8 and 
the  tliird  information  block  down. 

Finally,  tlie  label  pl/$  2+3  refers  to  page  1,  and  the  3*^^*  information 
block  alter  tlu'  label  $ 2.  p2/7,  8,  ')  refers  to  page  2 and  the  8*^', 

and  bl  ocks  . 


5.1,2  Description  of  Flow  Diagrams 

3. 1.2.1  Initialization 

a)  Central  Path  Angle;  pi /PROPAGATION  PRF.DICTION,  2,  3,  4: 

After  setup  of  the  registers  the  subroutine  TIIFTA  1 is  called  to  cal- 
culate the  earth  central  angle  between  the  submarine  and  the  station 
selected  by  the  argument. 


— ATAM 1 
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whcii’  I<| 

s’ 


Kfofontric  unit  voclor  directfd  to  the  submarine  estimated 
pos i t i on 

tlie  geocentric  unit  vector  directed  to  the  station 


and 


0-. 


TT 


b)  First  Magnetic  Parameter;  A^,  pl/5: 


Obtain  the  unit  vector,  AX,  normal  to  the  |iath  ol  transmission  from  tiu- 
station  to  tiu'  sulrnvirine  receiver 

s X r 

-AX  = — 

I S X it  J 


£li6: 

^1’  ^*2’  components  of  tlie  integration  position  P,  oriented  in  tlte 

geoientric  axis.  Tliis  vector  is  initialized  to  the  transmitter  position 
cooidinates,  tlien  moved  along  tlie  propagation  path  in  increments  ol  df/ 
raliiijis  as  the  integration  progresses. 

dt)  = 0.0)  radian  , 


nia- 

A^  is  The  first  magnetic  paranK'ter"  and  NP  is  a con.stant  geocentric 
unit  vectot  directed  toward  the  niagne tic  north  pole,  and  is  given  l>y: 


NP  = (0.9664,  0.0044864,  -0.2'.70S) 


Ca  leu  la  te : 


-(i\X  - N^) 
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) Sitisun.i  1 Index  IS;  |)2/2  throu^li  p3/() : 

'Hu-  st’asonal  ituli'x  IS  will  be  used  in  i^stimatin^  tbe  diuruni  uttfet 
later.  Theoretically  the  seasonal  year  can  be  sey;inenled  into  2A  2.16- 
week  piriods.  However,  t lie  currently  used  diurnal  t.ible  rerjuiri'S  on]-.- 
1 i ve  entries  . 


TABLK  I SKASONAI,  INDHX 


IS 

CALHNDAK  SPAN 

0 

JAN 

1 - 

.IAN 

1 5 

3 

JAN 

16  - 

.IAN 

31 

6 

FF.B 

1 - 

JUN 

30 

9 

Jin. 

1 - 

DKC 

15 

12 

DMC 

15  - 

DEC 

31 

d)  Sun  Vector;  S(I) , p4/$  2: 

Computation  oi  tbe  sun  vector  S(I)  = S(l),  S(2)  , S(3).  These  form 

a Kcecentric  unit  vector  representing  solar  information  required  to 
compute  the  diurnal  function;  i .e  . , a unit  vector  originating  at  the 
earth's  center,  pointing  at  the  sun. 

1)  To  start,  D , the  number  of  days  since  Jan  0.0,  of  1968  is  calcu- 
m 

lated,  to  bo  used  in  determining  the  sun  position,  and  from  that 
the  dayliglit  portion  of  the  propagation  path. 

This  can  be  accomplished  by  determining  the  number  of  days  in  prior 
years,  looking  up  in  Talile  2 the  number  of  days  in  the  months  pre- 
ceding the  current  one,  and  adding  these  to  the  number  of  days 
alreadv  past  in  tlio  current  month,  plus  the  fractional  part  of  the 
current  day. 
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TAtU.K  2 DAY  OK  THIS  YKAR 


(Him- lit  M 

Mont  h 

Ian  s!  1 

Keti  28  2 

Mar  31  3 

Apr  30  A 

May  31  5 

Jun  30  6 

Jul  31  7 

Aug  31  8 

Sep  30  <) 

Oct.  31  10 

Nov  30  1 1 


(,Ian  - 1 , Dec  = 


12)  Dm  (Days  in  prior  months,  since 
^in  1) 


Doc  31 


D = Number  of  days  since  Jan  0.0,  1068 


I)  = 365  (year  - 68)  + Dm  + (Number  of  days  in  current  month)  + MlAI’ 

t t 

(bast  two  digits) (From  table  above) 

^ (3600)  (hrs)  + (60)  (min)  + secs 

86400  (fraction  of  present  day) 

where  I.ICAP  = number  of  Feb  29  days  between  Jan  1,  1968  ami  present. 


2)  K,  the  fractional  part  of  D + 0.5 
K = (D  + 0 . 5) , modulo  1 


3)  I.^ , the  longitude  of  the  sun 


1,  = K,D  + K- 

s 12 


where  K^=  0.0172027914 


4.8649521 


represents  on  Jan  1,  1968  and  is  sublcct  to  update  to 


current  year. 
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'*')  M , ihf  orbital  anj;lo  of  tlio  earth  relative  to  the  perihelion  fpoinl 
s 

i;i  earth's  orbit  at  which  earth  is  closest  the  son). 

M = K I)  a K, 

S i a 

where  K = 0 .01  7201 'lb')') 

2 

K.  = 0.21030155 

i 

K,  rei>re.scnt  s on  .Jan  1,  1968  and  is  subject  to  update 
to  current  year. 

5)  First  compl  enKMit  ol  sun  vector 

S(l)  = K,  sill  L 
) s 

where  K,^  = 0.J978a36H  and  is  the  tilt  of  sun's  orbit. 

6)  1.,  sun's  apparent  longitude 

1.  = -2  7T  K - K.  sin  M - sin  (21,  ) 

6 s 7 s 

where  and  are  ellipticity  corrections 

0 / 

K,  = 0.0 134640 
0 

= 0.04127339 

7)  Second  complement  of  sun  vector 

2 ^ 

S(2)  = (l-S(l)  ) cos  h 

8)  Third  complement  of  sun  vector 

S(3)  = (l-Sd)'")  sin  L 


3. 1.2. 2 I’ath  Integration 

p4/$  55  through  i)13/l  constitutes  the  main  iteration  algori  thm  where  the 
Propagation  Prediction  terms  will  be  obtained  by  integration  in  i nc  remc'ntal 
steps  ol  0.01  radian. 

a)  Diurnal  Constant  Vector.  I’M;  p4/$  55  through  p5/$  15: 

Calculate  the  solar  zenith  angle,  Z,  from  which  the  diurnal  functions 
will  be  obtained  by  tlie  table  following. 

C().SZ  ■=  P • '.'3(1) 
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C.iUul.ili'  I h>  ilitiin.il  const, int  nsin^’.  I lu*  1 rctincnc  v , U',  .ind  the  season 

i ml*  X IS. 

TABl.K  1 DTI'KNAI,  CON.STANTS 


K (.'F,  IS)  IF  = frcqticncv  index 

IS  = season  index 


10.2 

11.3 

13.4 

0 

0.271 

0.241 

0.212 

3 

0 . 300 

0.248 

0.2  37 

() 

0 . 340 

0 . 304 

0.249 

9 

0.210 

0 . 1 Hb 

0.142 

12 

0.241 

0.214 

0.187 

NOTK : For  IS  = h ,ind  IS  = 9,  K is  at  present  constant.  It  is 

possLlile  that  this  would  not  he  so  in  a more  sophi  st  ic.U  cd 
version  ol  the  mode. 


tabu;  4 DIURNAL  FUNCTION 


TllKN  1 
IF 

! 

F (JF) 

FM  1 1 FM^ 

Cos  7.  < a 

1 

1 

1 1 

1 

0 

0 

a <Cos  7.  <b 

Cos  (Z-b)  + K(JF.IS) (a-CosZ) 

Cos  Z-b 

1 

a -Cos  Z 

0 

a-b 

a-b 

a-b 

b < Cos  7. 

K (.IF.  IS)  (1-Cos  Z) 

0 

l-Cos  Z 

Cos  Z-b 

1-b 

1 1 

L.  . 

1-b 

where  a = -0.15 
b - -0.04 

k - (IF, IS)  from  Table  3. 

and  Is  the  frequency  index,  1,  2,  and  3,  corresponding  to  10.2, 
13.f)  ad  11-1/3. 
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IS  is  t lu'  SI  .ison  index  detennintnl  in  t lie  description  ol  p2/2. 

FMj  , FM  j , i’M  ^ lie  Ihre.'  "sc.itter  dinrn.il''  Innclions  l liat  a|)ply  to  all 
tlirt'O  ! re<|iu’nc  ies  . 

F (.IF)  is  the  diurnal  liinclion,  and  is  calculated  for  each  trequencv. 


On  the  lirst  pass  through  the  diurnal  calculations  (above)  for  a ^ivcn 
station,  a value  l'(.JF)^  is  set  equal  to  FC.IF)  for  use  in  subsequent 

coiiiput  at  ions  . 


i)  Auroral  Compensation;  |i5/$  20  through  p7/$7  1 + 1" 

The  following  al gori t lira  wi 1 1 provide  compensation  for  auroral  effects 
when  ini  eg.rat  i ng  ever  a path  vector  near  the  pole.  For  values  of  con- 
stants refer  to  Table  5. 

The  lirst  calculation  is  the  magnetic  latitude  AA6 . 

A5  I P-iit'l  sin  O 

where  Nl*  is  a constant  vector  to  north  mag  pole 
AAh  - |aSIN  (A5)l 


pS/$  20-)  1 throuith  pb/7: 

The  auroral  constants  are  determined  by  the  following  iteration. 

If  AA^<AU(1)  1 = 1 

AU(1)<  M^<  AU(2)  1 = 2 

AU(2)  < AA^  <AU(3)  1 = 3 

AU(3)  < AA^  <AU(4)  1 = 4 

AU(4)<AA^  I = 1 

AAUR  = a (I)  * \A.  * AA^  + b (1)  * AA^  + c (I) 
bo  6 

BAUR  = Da  (I)  * AA^  * AA^  + Db  (I)  * AA^  + Dc  (1) 

b b b 
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TAHLK  5 CONSTANTS 


12  3 4 

1 1 1 


1 

1.02974 

1. 1344b 

1.16064 

1 

1.41371  1 

0 

310.04 

- 1805,55 

60.  30 

0 

- 641.39 

-f  4179.73 

-171.32 

0 

+ 331.71 

- 2414.90 

+121.65 

0 

+ 930,13 

- 5197,77 

- 93.80 

0 

-1905.09 

+11840. 17 

1 

+206.28  ; 

0 

+ 975.47 

- 6731.36 

-104. 13 

P .71: 

Finisli  fcilculat  ion  ol  aiiroi'.il  conslatils. 


B1  = -A2/ 


/77 


(AS)"  = sin  0 


B2 


I - 2 (Bl)  = cos  2Q 


B4  = z;  - AA^ 

i 0 


|0 


AA. 


= (B4)^  = ( 


AA^  = AA^ 


* Bl  = (^  - l^l  ) 


3 


s in  Q 


AA. 


. Bj  . ( i 

li- 


AA  «=  B4  * B2 
4 


) cos  20 
cos  20 
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FIGURE  2 CONDUCTIVITY  M/\i> 
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' ^ riu’  Ci>iuliii- 1 i \' i 1 ItkIi'x,  TC;  pH/1  ilii'ou)',li  i)l()/$ 

This  .11^01  itlun  takes  tin  position  vt-ctor  1’  and  snarclu's  the  conduct  ivitv 
m.ip  to  tind  tile  conduenvLty  index  J C of  the  rff>i  on  of  tliP  cartli's  siir- 
t.u  c I' r ns  piMid  i nj;  to  i’. 

Till-  ci'ndiu- 1 i vi  t V map  (Figure  2)  has  t lu'  following  lorm: 

F.'f  latitude-  north  ot  the-  euirth  is  hroke-n  eip  into  hoxe'S  e-ach  s” 

h 

F >)!  latitueles  soutli  e)l  ^ ■)  S,  the  e'artli  is  hrejkcn  up  into  lioxos  5^  in 
latitueh  anel  10  in  1 e'llgi  t Uele’  . 

There-  are-  1 '>  coiiehu- 1 i vi  t v le'vels  ranging  from  16  (sea  wate'r)  to  1 
( pe  rmaf  ros  t ) . 

The'  ceuuhic- 1 i vi  t V e>f  "ach  hox  is  the  average  coriduc  t i vi  t v of  the  e-arth's 
surface'  within  the  hox,  weighted  accoreling  to  the  areas  of  each  con- 
ductivit'.-.  The'  Wi'st  ing.house-  Conductivity  Majj  was  useel  . 

The  three-  pages  which  eh'serihe  the  conductivity  search  are  straight- 
ferrwarei  euiel  eie)  llejt  re-ejuire  e' 1 ahorat  ion  . 

el)  IJiurnal  Funetie)ii,  K;  pi  1/3,  4: 

From  Taliie  T uneh  r de'seription  for  p4/$56  the  valeie  for  K(,IF,IS)  is 
ohtainiel.  The-n 

F Rl^  + K (.IF,  IS)  m., 

whe-re  IF  is  use-el  .is  index  to  distinguish  it  from  the  diurnal  function. 

If  this  is  first  tinv  through,  then  go  to  $A1 ; otherwise  continue, 
e)  Correction  Facteirs,  ^'21'’  P''/'’- 


C^j  is  the-  corre-ction  factor  lor  reconverted  first  mode  of  propagation. 
Cll^.JF)  = C^jf.IF)  + Kjj(.IF)  (F(.!F)-F(.JF)^) 

ti^j  is  the  correction  factor  for  higher  propeigation  modes. 

C^|(.IF)  = C^jf.IF)  o^(JP)  K(.IF)  + K^jdF)  (F(.IF) -F(.IF)  ^^) 

He  re  se-t  Ff.IF)  - F(.IF) 

o 

wh.ere  K^j  and  j are-  constants,  defined  for  each  frequency,  and  repre- 
sent the-  excitation  of  2nel  mode  and  re-e-xc  i I .U  ion  of  Is  I mode,  o is  the 
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ri' 

1 .it  i V' 

night  t i mi- 

at  t I'nuat  ion  ol 

2nd 

moili' 

and  i s a 

1 unc  l 

i on  o f 

fro 

(JUflU  • . 

f> 

(1) 

U 

0 

.91 

Kii(l) 

= 6 

.7 

^21 

M) 

= 25. 

0 

o 

(2) 

o 

0 , 

4)2  5 

Kii(2) 

= )i 

.7 

^'21 

(2) 

= 25. 

0 

o 

(■}) 

o 

0 , 

. ‘)4 

K^l(3) 

= h 

.7 

^21 

(3) 

= 25. 

0 

i')  Diurn.il  I'liiKtion  AviTa:,',».' , CC\  pll/$  42: 
('.(UJI-)  = (:G(.in  ♦ KMK) 

1 ) I nt(.',r,r.U  i on  ; i)l2/l,2  : 

h 

g - N (;AY(.I1M)Mj  + CAYl  f.IF.IC) 
1-2 


AUK)  Q + AAUR* 

()  = '^  I)c:AY(.IK,I)  t-  DCAYl  (JIMC) 

1=2 


lU.Jl-')  = ')  c;,^  f 15AUR* 

A()„  = C,  C,  C.  RATfJK)I>f 

2 14  5 1 


StTM  A (>2  = ' AO., 


AO-j  - (A(.IK)  + F(JF)  B(.JF))  R.\T(JF) 


SOM  AO.j  = S^IM^0,J  +^0^ 


whirc  t;AY , CAYl,  DCAY , DCAYl  can  be  found  in  tables  6 and  7. 

Cj  is  tite  ellipsoidal  correction  constant 
= -0.336  X in"^ 

C,  is  the  chart  conversion  factor 
4 

= 0.6974 

is  the  nominal  10.2  wavelengths/radian 
= 21/ .021 . 

*NoLe  that  in  the  1 low  diagrams  the  terms  AAUR,  llAUR  were  erroneously  excluded. 
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-32  .0 
-27.0 


-4U  . 


- 8.0 

3() . 3 

32.7 
25.3 
19.1 

14.8 

12.0 

9.2 

7.5 

6.2 

5.2 

4.5 

3.8 

3.9 


-35.0 

-23.0 

27.5 

44 . 0 

38.5 

33.2 

30 . 0 

27.6 

25.6 

23.8 
22  9 

22.2 
21.5 
21  .0 

20.8 


IC^  10.2 

11.3 

13.6 

1 

0.0 

30.0 

60.0 

2 

3.8 

28.2 

55.0 

3 

8.3 

21  .2 

38.5 

4 

5.5 

5.5 

12.0 

5 

22.3 

19.4 

15.2 

6 

26.3 

23.3 

19.2 

7 

30 . 2 

27  .0 

23.3 
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TARI.K  7 (c  out  i iiui'(i) 


it^ 

10.2 

11.3 

13.b 

8 

32  .8 

20.5 

25.2 

9 

3 3 .8 

30.0 

25.2 

10 

14.2 

30.2 

25.2 

1 1 

3 3.b 

10 .0 

25.2 

12 

3 3.2 

20.4 

24.4 

13 

32.0 

28.0 

23.7 

14 

32  .b 

28.8 

23.7 

15 

32.1 

28.5 

23.7 

lb 

31.8 

-ft 

27  . b 

23.7 

Multiply  all  by  10 

ft)  IntcKration  ComplPtc  Test;  pl2/4  tfirousb  plJ/I; 

The  intcj’ration  counter  sums  the  number  oi  dO  increments  used  where  dO  = 
0.01  radian. 

Ti:ST  = 

If  ] flTsl  j > (d(i)^ 
tlien  Increimmt  1’ 

P = (;\X  X ■?)  • Sin  (do)  + T>  • cos  (dO) 

Return  to  Paragraph  p4/$  55. 

If  I Test  I ^ (dO)^ 

Then  continue. 

3. 1.2. '3  Generation  of  Output  Data 
a)  Averat;inK  Constant.  AVG;  pl3/$  ()(H1: 

Corrects  lor  inc  ri'mt'tit  al  intenration 
AVG  = Ol/n.V, 
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h) 


O 


cl) 


win  r<  II  is  iiuniHT  ol  i ncrtT*' lU  s uscil  and  AO  = 0.01  radian. 
li.isi'  Station  (aiiti.il  An>;l>  : 

(.ill  sill' riMit  i 111  TllKTAl  to  compiiti-  Oj(haso) 

.Sui'iiii.il  ion  Avi  rai',1  ; pi  4 / 1 : 

ASl-M  0.,  AVC.  (SlT-no,,) 

A.srM  0^  = AVC  (srM\o^) 

(1  rt'ijul'ncv  di'pendoncc  implied) 
rroiKit;at  ion  Oorrect  ion;  Oj^ 


0.j(.]F)  = O.^C.JK)  -t  I'.AV  (.II- , 


1) 


ViGC(.lF)  DCAY  (.IF,  1) 
f'^(.JF) 


e)  Predicted  Variance;  „ .;  pl.4/3: 

•'  Pill 


Comi'ine  tlie  el  feet  s ol  second  mode  and  reconverted  first 
estimate-  of  the  variance  on  the  iiredicted  phase. 


2 

0 

ppi 


(c^i  -t 


f .04)' 


to  compote  an 


f)  .Spheroidal  Correction; 


(“^(.IF)  = ASIT-U)^ 

'3. 1.2. 4 Tracking  Filter  Update 

a)  Predicted  Tracking  Filter  Measurement,  6p(.lF)  : 

pi 7 

fi^,.(.JF)  lor  station  i is  obtained  by  subroutine  T.IF.TAC . It  is  the 
predicted  phase  lor  station  i. 

, (.IF)  for  base-station  k is  similarly  calculated. 

Ck 
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-„(jn  = ■■  . (M-)  - - , (IK) 
I’  ^ t 1 c-  k 


p 1 5 / $ 't  S t-  1 : 


The  iiueger  Lane  counl  is  adjusted  in  the  tracking  t liter  for  station  ik 
on  Irequency  .IK. 

P'W$  nS  + 2: 

It  more  lreciucnci.es  lo  process,  increment  IK  and  return  to  pn/$  hi, 
olhe  rwise  exit  . 


'3.2  Fl-tiW  DIACdUM.S 

Tin  I’rop.igat  i on  Prediction  Subprogram  flow  diagrams  are  presented  on  the 
following  pages. 
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PROPAGATION  PRCOICTION 

This  PROCpRAM  will  compute  the  predicted  spheroidal  at*  propagation 
CORRECTIONS  TO  the  NOMINAL  predicted  PHASE  MEASURCTCNTS  FROM  A GIVEN 
station  it  also  COMPUTES  AN  ASSOCIATED  VARIANCE.  TTC  Tt^TEE  PARA- 
TETERS  ARE  COMPUTED  TOR  ALL  J EREOUENCIES  AND  SAVED  IN  MEMORY  FOR  EACH 
STATION  THE  total  LANE  COUNT  FROM  TmE  CORRECTED  PREDICTED  PHASE 
REPLACES  ANY  lane  counts  CONTAINED  IN  TME  ECaSUREO  PHASE  IN  T*« 
tracking  eilter  ttc  Staticjn  Nur«ER  IS  tec  OH-r  AROOrCWT. 


I E>ROPAOATJON  PRCDtCTJ 

_ _ i 

CREATE  NS  WORDS  OF  SCRATCH 
storage  in  the  R|5  POSMOOAPE 
stack.  ZERO  ALL  MS  WORDS 


• NOTE  - THE  calling  PROGRAM 

• WILL  HAVE  TO  DELETE 

• THE  MS  WORDS  FROM 

• THE  R15  STACK 


0001 
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pAce  0002 


DOOM 
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Page  0003 


This  page  left  blank  intentionally 
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p*GC  000s 


■—n- 

COMT-UTC  SON  POSITION  VECTOR 
S3  - 1 1 I-ISI •••ai**  SI'COSILT 

sa  • 1 1 1-isi •••ai** .5i*siNiLi 

L - -aiPI IE  ♦ K651NIMSI 
- KTSlNiacSI 

SI  - KSSINILS* 

MS  - K3U  ♦ KS 

LS  • KID  ♦ « 

E - fractional  part  or  0 ♦ .s 
0 - nu>w;r  or  days  since  1/1/72 
SM)  scaled  BI 


• START  integration 

• FIRST  COrr^JTE  TERMS  THAT 

• ARE  ITCtPENDANT  or  FREQUENCY 
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Fn^  • (A  - COS  Zi/iA  - ai 
rM3  - 0 


• COMPUTE  AUPOPAL  TERMS  AAIII 

• I • a TO  6.  AAUR  Arc  BAUR 

1 

«0 

’ COM^TE  AA6 

, A5  - ABSOLUTE  IP  . magnetic  POLE 
VECTOR!  scaled  B1 

AAG  • ARCSIN  AS  AT  PIB-I 

i IAA6  IS  the  magnetic  L at  I 




SET  UP  TO  COMPUTE  AAUR  AND  BAUR 

WHICH  ARE  A Function  or  the  mao 
NETIC  latitude 

SET  AAUR  - 0,  BAUR  - 0.  1-0 


I TEST  AA6  against  the  SET  OT 
L I > TO  DETERMINE  I 


«72 

FORM  AA6  - AU  LIMIT  table  III 

ooos 
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PAGC 


1 




COMPOTE  AA^r 

AA3.  AAS  and  AAS 

AA9 

• BM»«3 

AT  Pl*»3e-1 

AAM 

• Bs»B2 

AT  PI 

AA3 

- AA5*B3 

AT  PI**38-2 

AAa 

- AA5*B1 

AT  Pl*«38-2 

- Plr^  - 

AA8 

- I - e(Bn»»2 

Bl  ■ 

-A2/  1 1 

A5**a  •••.s 

SAve  AA^  AA3  AA4  AA3  AA6 
AAue  Af«  BAUR 


0007 
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OtTERHirC  COKWCTIVITy  LEVCL 

The  rOLLOMlNO  section  or  the  propagation  prediction  program  uses  Trc 
conductivity  map  to  determine  TK  CCPCmCTIVl  Ty  level  at  the  present 
point  on  the  path  or  integration  OCriNEO  BY  TK  VECTOR  P.  TOR  PRO- 
GRAMMING ECriCIENCY  TK  EARTm  HAS  BEEN  DIVIDED  INTO  BOXES  EVERY  rivt 
DEGREES  or  LATITUDE  AND  LONGITUDE  EXCEPT  THAT  BELOW  SOUTH  SS  LATITUDE 
The  longitude  SIDES  ARE  10  DEGREES  IN  LENGTH,  TrC  MAP  HAS  BEEN  DIVI- 
DED INTO  J6  PARTS  (1  rOR  each  INCREMENT  OT  LATITUDE  I AND  A TABLE  OT 
ADDRESSES  POINTS  TO  THE  PART  OT  Tf<  MAP  THAT  IS  TO  BE  USED  TOR  A 
GIVEN  latitude  the  proper  ENTRY  IS  OETERMirCO  BY  COMPUTING  THE 
NUMBER  or  5 DEG  LATITUDE  INCRETCNTS  ABOVE  THE  SOUTH  POLE.  TiC  TIRST 
Entry  is  rOR  90S  to  BSS  ANO  the  last  entry  is  for  BSN  TO  90N.  FOR  A 
GIVEN  latitude  THE  H BIT  CONDUCTIVITY  LEVELS  FOR  THE  LONGITUDE  INCRE- 
MENTS ARE  STORED  IN  CONSECUTIVE  rCXAOECIMAL  CHARACTERS  STARTING  AT  20 
DEG  EAST  longitude  AND  CONTINUING  ALL  THE  WAY  AROUND  THE  EARTH  FOR 
increasing  east  longitude.  the  CORRECT  CONDUCTIVITY  LEVEL  IS  THEN 
OBTAINED  BY  COMPUTING  THE  NUMBER  OF  LONGITUDE  INCREMENTS  AT®  COUNTING 
THWCXXJH  The  TABLE  UNTIL  THE  CORRECT  ITEM  IS  REACKO. 

EACH  16  BIT  computer  WORD  CONTAINS  THE  CONDUCTIVITY  LEVEL  FOR 
FOUR  increments  or  longitude.  TK  most  significant  H bits  or  THE 
first  word  or  THE  MAP  FOR  A GIVEN  LATITUDE  CONTAIN  THE  LEVEL  BETkCEN 
20  AND  2S  DEG  EAST  LONG.  THE  LEAST  SIGNIFICANT  M BITS  ARE  FOR  LONG 
E35  TO  EHO.  the  MAP  THEN  CONTINUES  IN  THE  NEXT  HIGHER  MEMORY  LOCA- 
TION FOR  the  next  longitude  INCREMENT.  THE  CONDUCTIVITY  LEVEL  IS  A 
number  between  0 AND  IS  THAT  IS  USED  TO  INDEX  A TABLE. 

BECAUSE  MOST  or  THE  EARTH  IS  MATER  AND  THEREFORE  HAS  TK  SAK  CON- 
DUCTIVITY LEVEL  A COMPACTING  TECHNIQUE  IS  USED  THAT  REDUCES  TK  SIZE 
OF  THE  TABLE  BY  APPROXIMATELY  HALF.  WKKVER  TK  LEVEL  FOR  WATER  (131 
IS  ENCOUNTERED  F(3R  A LONGITUDE  TK  KXT  KXADECIMAL  CHARACTER  IS  A 
COUNT  OF  TK  NUMBER  OF  CONSECUTIVE  BOXES  OF  WATER  FOR  THAT  LONGITUDE 
UP  TO  A MAXIMUM  or  15. 
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PACC  0009 


cmangc  Lonoiruoe  to  to  otCRtt 

INCRt^CNTS 

931 

SET  UP  FOR  S 

READ  ADDRESS  OF 
LATITUOt 

ESTABc ISM  S COUN 

number  of  CmaRAC 
SET  THE  COUNTER 

A NEW  UORO  IS  NE 

EARCM 

HAP  FOR  THIS 

TER  TO  COUNT  T»C 
TCRS  IN  A MORO 

TO  ITCICATE  that 
COCO  ireCOIATELT 

• start  search 
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PACC  00 1 0 


CLtAB  PCOlSTtP  to  MAKE  ROOM  ro* 
IMC  COUNT  Of  CONSECUTIVE  LEVELS 

or  mater 


SHIFT  WORD  left  ONE  DIGIT  TO 
PLACE  TFC  COUNT  OF  CONSECUTIVE 
LEVELS  or  hater  in  A SEPARATE 
REGISTER 


FORM  longitude  I 

- ref>eat  count 

IREF>EAT  COUNT  - 
CUTIVE  BOXES  - 1 

NCREMENTS 

NO.  OF  CONSE- 

) 

RESTORE  CONO  LEVEL  OF  WATER  THAT 
WAS  LOST  IN  GETTING  REPEAT  COUNT 

S36 

DECREMENT  THE  NO 

INCREMENTS  BY  1 

. OF  LONGITUDE 

IS  THIS  TTC 
:;:iEvEL  FOR  the  given  LONGITUDE? 

( INC  - 0 OR  ->  ^ 


NO 


0009 


rzs 

_ _ 3 

1 SAVE  the  conductivity  level  no 

that  part  or  Ttc  INTEGRATION  THAT 
IS  itcepenoant  of  frequency  IS 
complete,  non  start  terms  that 

ARE  A function  OF  FREQUENCY. 
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• TEST  FOR  END  OF  INTEGRATION 


_J 

DECREMENT  INTEGRATION  COUNTER 
BT  .01  RADIAN 

.01  RADIAN  IS  the  integration 
INCREMENT 

integration  counter  is  THE  PATH 

length  remaining  for  integration 
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PAGE  00 1 1 

rotate  p Through  t^c  integration 

INCREMENT  OA  - .01  RADIAN 

PINEH)  - PCOSD*  ♦ lAX  X PISITCA 

1 

SM 

DOOM 


INTCWATJON  Ct3r»n.CTE 

NOTE  - the  integration  mas 
actually  perfoko 

OVER  A RANGE  OF  N T|KS 
THE  INTEGRATION  INCRE- 
MENT or  .01  RADIANS. 
thetai  is  less  than 
OIN  LUICM  IS  LESS 
THAN  thetai  ♦ .01 
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e*CC  00 1 >4 


• propagation  prediction  is  now 

• COMPLETE  FOR  THIS  FREOUCNCV 

• PREPARE  TO  UPDATE  LATC  COUNT 

• IN  TIC  tracking  FILTER  PHASE 

• estimates 
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T»C  ARCUrCNT  IS  UAOISTURKO 
^S  MOnOS  HAVC  BCCN  AOOCO  TO  T»< 
RIS  POSH  OOM4STACK 
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Kf  COMl’rTKK  SinU’R(H:i<AM  KNV''!  RON'MI'.NT 


i . i . 1 l’i'ii|i.i)’,.it  ion  I’rt'd  i c t i on  T.iblcs 

al  M.i);ni  l ic  I’oli  : Tliis  table  contains  the  t h roe -e  1 emi'iit  vector  that  tIcfiru'S 

t tu’  position  ol  the  Nortli  Magnetic  Pole.  it  is  defined  in  detail  in  the 

listing.  The  official  nam.'  of  this  table  is  NORTH  MACNKTIC  POLK  VKCTOR. 

b)  Auroral  Oompensat  ion  Limits:  Thi.s  table  is  used  to  test  the  Magnetic 
Latitude  lor  selection  ol  the  appropriate  auroral  constants.  The  first 
eiit  rv  i.s  the  lowest  limit  (AU(1)),  the  second  entry  is  AU('2) -AUf  1 ) , the 
tliird  eiit  i".'  is  AU(3)-AU(2)  and  the  last  entry  is  AU(4)-AU(3).  .See  the 
listing  for  more  detail.  The  official  name  of  this  table  is  AUR  LIMITS. 

cl  Oonduc t i vi t V Table:  This  table  defines  the  conductivity  level  for  any 

point  on  the  earth.  It  is  defined  in  detail  in  the  fleiw  charts  anei  in 

the-  listing.  rile'  official  name  of  this  table  is  OONI)  TABLK. 

el)  Diurnal  Constants:  This  table  contains  the  diurnal  constants  for  the 
calcul.it  ion  of  the  diurnal  function  F.  For  convenience,  the  five 
different  entries  in  the  original  table  for  the  24  values  of  the  season 
inele-x  h.ive-  bee-ii  renumbered  0 to  4.  This  table  is  defined  in  eletail  in 
the-  listing.  'Hie  official  name  of  this  table  is  DIURMAL  CONSTANT  TAIU.K . 

e)  Second  Moeb-  Kxcitation  Factor:  This  table  is  used  in  the  computation  of 
PROP  C.  It  is  lie-line'd  in  detail  in  the  listing.  Tlie  official  name  of 
this  t .lb  1 e is  K2 1 . 

i)  Second  Mode  Relative-  Night  Attenuation  Factor:  Tliis  table  is  used  in 
the  computation  of  PROP  C.  It  is  defined  in  detail  in  the  listing. 

The  official  name  of  this  table  is  ALPHA  0. 

g)  Model  Constants:  This  table  contains  five  entries  for  CAY(2)  to  CAY(b) 
used  in  the  computation  of  A.  Kach  e-ntry  has  been  modified  by  C4  and  C5. 
It  is  defined  in  detail  in  the  listing.  Tlie  official  name  of  this  table 
is  CAY  TABLK. 

h)  Model  Constants:  This  table  contains  five  entries  for  each  frequency 
for  I)CAY(2)  to  DCAY(6).  Kach  entry  has  been  modified  by  C4  and  C5 . 

It  is  defined  in  detail  in  the  listing.  The  official  name  of  this  table 
is  DCAY  TABLK, 

i)  Conductivity  Constants:  This  table  contains  lb  entries  for  each  frequenev 
for  CAYK'I')  to  CAY1(16).  Kach  entry  has  been  modified  by  C4  and  CS. 

It  is  defined  in  detail  in  the  listing.  The  official  name  of  tills  table 
is  CAYl  TABl.K. 
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i ")  c:oiuhi(  tivity  Const.tiUs:  This  l.ihli'  cniif.iiiis  Ih  i-ntriis  lor  oach  freqiioncv 
lor  IXlAYlil)  lo  I)('AY1(1<>)  . Kacli  ciu  rv  lias  hi’cn  rnmlilii’d  by  ' 4 and  fY) . 

It  is  dolini'd  in  de  tail  in  thf  listiny;.  Tlif  official  name  of  this  table 
is  IX'AYl  TAHhi;. 

k)  Frequency  Conversion  Factor:  This  table  contains  the  factors  that 
convert  10.2  d.it  a into  1)  . b or  1 1 . i . It  is  defined  in  detail  in  the 
listinj;.  Tliis  official  narx'  of  this  table  is  RAT. 

l)  F.llipsoidal  Oorrect  ion:  This  table  contains  Cl  modified  by  C4 , O')  and 
the  iretiuencv  conversion  factors.  It  is  defined  in  detail  in  the  lisl- 
i ny, . The  official  name  ol  tliis  table  is  Cl  TABLK . 

m)  Model  Constants:  This  table  contains  (2\Y(1)  modified  by  C4 , C5  and 

the  frequency  conversion  factor.  It  is  defined  in  detail  in  the  listiny. 
The  official  natiK’  of  this  table  is  CAY  ONK . 

n)  Model  Constants:  This  table  contains  DCAY(l)  modified  by  C4 , C'l  and 

the  frequency  conversion  factor.  It  is  defined  in  detail  in  the  listiny. 
The  official  name  ol  the  table  is  DCAY  OKK . 

o)  Auroral  Constants:  This  table  contains  the  three  sets  of  auroral 
constants  modified  by  C4  and  C5.  It  is  defined  in  detail  in  the  listiny. 
Tlie  official  name  of  this  table  is  AUR  DATA. 

p)  .Station  Locations:  Tliis  table  contains  the  locations  of  all  exist  iny 
OMKGA  transmitting  stations.  Kach  location  is  specified  in  a tliree- 
element  grocentric  position  vector.  The  first  entry  is  for  station  A. 

The  table  is  di'lined  in  detail  in  the  listing.  The  official  name  of 
this  table  is  STATION  VKCTOR  TABLF . 


'3.3.2  Propagation  Prediction  Temporary  Storage 

All  temporary  storage  used  by  the  Propagation  Prediction  Program  is  in 
the  R15  pushdown  stack,  except  for  seven  words  associated  with  the  Auroral 
corrections  called  AUR  SAVK  which  are  placed  iti  the  Kalman  program  storage 
for  new  state  vectors  during  measurement  update. 
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. 5 . 3 Ri’qiiirod  Syslfm  Library  Subroutines 


SI’HROITINK 

Ki.ow  DiACHAM  rkfkri;n(:k 

SUBPROGRAM  DKSIGN 
DOCUMI3NT  (by  Volume  Number) 

AS  IN 

l>')/$  20 

CROSS 

pl/i 

All  Subroutines  found  in 

I’l/^ 

Volume  XII 

, DOT 

pl/3 

SIN -Cos 

p4/l 

SQRT 

pI/5  p7/$  71 

1 

p4/l 

1 TllKTAl 

1 

pl/3 

j THKTAC 

pl4/() 

. 

1 

pl5/2 
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